arXiv:1506.02590v2 [physics.class-ph] 20 Sep 2015 


What’s the Use of Bound Charge? 
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Bound charge is a useful construct for calculating the electrostatic field of polarized material, and 
it represents a perfectly genuine accumulation of charge. But is such a material in every respect 
equivalent to a particular configuration of bound charge? The answer is no, and the same goes for 
bound current and (in the time-dependent case) polarization current. 


I. INTRODUCTION 

In introductory electrostatics we learn that the elec¬ 
tric field of a polarized object (polarization P = electric 
dipole moment per unit volume) is equivalent to the field 
produced by surface and volume “bound” chargei^^ 

CTf, = P • n, pb = -V • P (1) 

where n is a unit vector perpendicular to the surface 
(pointing outward). This is easy to understand: po¬ 
larization results in perfectly genuine accumulations of 
charge,^ differing from “free” charge only in the sense 
that each electron is attached to a particular atom.^ 

But is polarized material in every respect equivalent to 
such a distribution of bound charge? For example, is the 
electric force on a polarized object the same as it would 
be on pb and Cb? What about the torquef And how about 
the force and torque densities within the material? 

The very notion of density (whether of force, torque, 
energy, or even mass, charge, and dipole moment) can be 
problematic. After all, matter is composed of atoms, and 
on a microscopic scale these quantities fluctuate wildly in 
position and time. We mean, however, their macroscopic 
averages over regions large enough to contain enormous 
numbers of atoms and yet small compared to the relevant 
dimensions of the object. 

We will confine our attention to classical macroscopic 
electromagnetic forces. Of course, the atoms in a solid 
or liquid are subject to all sorts of “mechanical” forces 
(which may themselves be electromagnetic on a micro¬ 
scopic level), and they are governed by the laws of quan¬ 
tum mechanics. But in this paper our purpose is to ex¬ 
plore the role of bound charge, and to this end we adopt 
a radically simplified model: We imagine a continuum 
of ideal (neutral) point dipoles, described by a specified 
function P. How this polarization came to be, and what 
“mechanical” forces sustain it, we do not inquire. (Imag¬ 
ine, in the static case, that they are simply glued in posi¬ 
tion.) We are interested only in the electrical forces ex¬ 
erted on these dipoles by the (macroscopic) field E (the 
total field, attributable both to the dipoles themselves 
and to any external sources). 

If the question is “How does a particular deformable 
medium respond to externally applied fields?” then one 
requires detailed information about the structure of the 


material, its elastic and dielectric properties, the pres¬ 
sure, the temperature, and so on.® Our question is much 
simpler: “For a stipulated polarization, what is the elec¬ 
tromagnetic force density, and in particular can it be 
calculated by replacing P with the associated bound 
charge?” We take it to be the “correct” force density 
(as distinct from the force density associated with bound 
charge), but remember that it does not include the “me¬ 
chanical” stresses that would also be present in any real 
material. 

In Section H we rehearse the standard derivation of the 
electrostatic potential of a polarized object, in terms of 
the bound charge. We then apply the same reasoning to 
the force and torque on the object. In Section HI we do 
the same for static magnetization, and in Section IV we 
generalize to time-dependent configurations. In Section 
V we compare our results with the Einstein-Laub force 
formula, and in Section VI we draw some lessons and 
conclusions. 


II. FIELDS AND FORCES FOR POLARIZED 
MATTER 


Let’s review how bound charge is first introduced: The 
potential of an ideal dipole p is 


V(r) 


1 p • 

47reo ^ ^ 


( 2 ) 


(where = r — r' is the vector from p, at r', to the field 
point r). The potential of an object with polarization P 
is therefor^ 


V(r) 


1 

47reo 


P(r') • i 
^ 2 




The standard integration by parts, using 


V' 




( 3 ) 

( 4 ) 


turns this into 


V = 


1 

47reo 



V'-P) 

n- 



P ft 


n- 



( 5 ) 
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and we conclude that the potential of a polarized object 
is the same as that produced by the charge distribution 
Pb and (Tb- 

Next we ask “What is the force on a piece of polarized 
material, in an electrostatic field E?” The force on an 
ideal dipole p in an external field E i^ 

F = (p. V)E, (6) 

so the force on a chunk of polarized material is 


What about the torque on a polarized object, in a static 
electric field? The torque on an individual dipole i^ 

N = (pxE) + [rx (p-V)E], (12) 

where r is the vector to p from whatever point we choose 
to calculate torques about. The net torque on a polarized 
object, then, is 

N = y (P X E) + y [r X (P • V)E] d^r. (13) 


F= / (P-V)Ed3r. (7) 


As before, we integrate by parts; the Ah component is 



Thus the total force on the object is the same as it would 
be on pb and CTb- 

However, the force densities are not the same. Equa¬ 
tion says the force per unit volume i^ 

f=(P-V)E, (10) 

whereas Eq. suggests 


As always, we integrate by parts: 

A^i = y Sijk [PjEk + rjPi{ViEk)] d^r 

= I e.jk[P,E, + Viir,m,) 

- (yirj)PiEk - rj{\7iPi)Ek)]d^r 
= y {V • [P(r X E),] - (V • P)(r X E) J dh, (14) 

(summation over repeated indices implied; = Sij), 


N = (r X E)(P • n)da - y (V • P)(r x E,)d^r 

[r X (pbE)] d^r. (15) 


= J) [r X (cTbE)] da - 


Again, the total torque on the object is the same as it 
would be on the bound charges. 

However, Eq. 13 indicates that the torque density in 
the material is 


n = (P X E) -I- (r X f). 


whereas Eq. 15 says it is 

nb = r X PbE = (r X fb). 


(16) 

(17) 


fb = PbE = -(V-P)E. 


( 11 ) 


These two expressions are certainly not equivalent. Imag¬ 
ine, for example,^ a “bar electret” (a cylinder uniformly 
polarized along its axis); pb = 0 (the only bound charge 
resides on the two ends), so fb = 0 , but if the field is 
nonuniform f 7 ^ 0. Bound charge incorrectly distributes 
the force, even though it gets the total force right. 

This raises a surprisingly delicate question: What do 
we mean by the “force density” inside the medium? Pre¬ 
sumably we should (in the mind’s eye) isolate an infinites¬ 
imal piece, of volume u, determine the force on it, and 
divided by v. But this little piece carries surface bound 
charge in addition to its volume bound charge, and it’s 
easy to see (reading Eqs. ^ in reverse) that the total 
force is precisely (v. Of course, in the bulk material the 
surface charge on v is canceled by that on the adjacent 
inner surface of the surrounding medium—there is no net 
“surface” charge within the substance—but if we’re inter¬ 
ested in the force on v alone, its surface charge must not 
be ignored. The force density fb (Eq. is incomplete, 
because it does not include this contribution. 


These expressions are not equivalent. Eor example, if P 
and E are uniform, nb = 0, whereas n = (P x E) —and 
surely there is a torque on the dipoles. Once again, treat¬ 
ing the medium as a configuration of bound charges gets 
the total right, but incorrectly assigns its distribution, 
because it ignores the role of “internal” surface bound 
charge. 


III. MAGNETIZED MATTER 


Now consider the magnetostatic analog: a chunk of 
magnetized material (M = magnetic dipole moment per 
unit volume). The vector potential of an ideal dipole m 
is 


A = 


po m X -4. 


Itt ^ ^ ’ 

so the potential of the magnetized object is 




M(r') X ^ 3 

>,2 1 , 


(18) 


( 19 ) 
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and integration by parts (again using Eq. yields 
(V' X M) , / (M X ft) 


A=^ 

47r 


d^r' + 


Uv 


da' 


( 20 ) 


The two terms are identical to the potentials of (bound) 
volume and surface currents P 


Jfe = V X M, Kfa = M X ft. (21) 

Once again, these are perfectly genuine currents, differing 
from free currents only in the sense that they are the 
collective effect of many tiny current loops—as in a relay 
race, no particular electron makes the entire trip. 

But is magnetized material in every respect equivalent 
to the currents Jf, and K;,? For example, are the forces 
on them the same? The force on a magnetic dipole m 
iPl 


F = m X (V X B) + (m-V)B, (22) 

so the force on a chunk of magnetized material is 

F = y [M X (V X B) + (M • V)B] d^r. (23) 

From the vector identity 

V(M • B) = M X (V X B) + B X (V X M) 

+ (M • V)B + (B • V)M (24) 


it follows that 


F = 


V(M • B) - B X (V X M) - (B • V)M 


= / (Jf, X B)d^r + G, 
Jv 

where 


d^r 


(25) 


G,= [ [V,(M,B,) - Bj(VjM,)] d^r. (26) 

Jv 

The second term in the integrand is Vj{BjMi) — 
= VjiBjM,), scP 

G, = / [VdMjBj) - V,(M,B,)] d^r 
Jv 

= y \hi{MjBj) — nj{MiBj)\ da 
= ® [(M X ft) X B]^ da 


Js 

= y [Kf, X B] . da. 


(27) 


Thus 


F = / (Jfe X B) d^r + (i) (Kfa X B) da, (28) 


However, the force densities inside the medium are dif¬ 
ferent: the force per unit volume on Jb would be 


fb = (V X M) X B, 


whereas the force density (from Eq. 231 is 

f = M X (V X B) -t (M- V)B. 


(29) 

(30) 


The torque on a magnetic dipole m in a magnetostatic 
field B itP 


N = (m X B)-br X [m X (V X B) + (m-V)B]; (31) 
the torque on a magnetized object is therefore 
N = y |(M X B) 

-br X [M X (V X B) + (M- V)B]}d^r. (32) 
Using the identity 


^pqr^pst — ^qs^rt ^qt^rs^ 


(33) 


A, - y 

+ Mi(ViHfe)]}d4 

= y Gjk [MjBk + rjM;(VfcH;)] d^r (34) 

= y e^jk [MjBk + \7k{rjMiBi) - rjBi(ykMi)\ d^r. 

Subtracting and adding 

Vi{rjMkBi) = {Virj)MkBi + rjiyiMk)Bi -b rjMkiViBi) 

= MkB,+rjBi{ViMk) (35) 

to the expression in square brackets (last line of Eq. , 
we find 

N. = J,,,{lVUr,Am)-Vdr,MM)] 

+ r,Bi[{ViMk)-iVkMi)]}d^r. (36) 

We are now set up to integrate by parts, using 

/ (VfcQ)d^r= f Qfikda, (37) 

Jv Js 

where the function Q may carry one or more indices. 
Thus 


Ni = f e^ijkTjBi [[Ml hk) - [Mk ft/)] da 


+ / e^ferjH/[(V/Mfe) - (VfcM/)] d^r, (38) 
Jv 


and so 


and the total force on the object is indeed the same as it 
would be for the bound current distributions. 



[r X (Kb X B)] da+ [ [r x (Jb x B)] d4. (39) 

Jv 
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Once again, the bound currents get the total torque right, 
but whereas the torque density (from Eq. 32) is 


n = (M X B) + (r X f), 
the bound currents suggest (Eq. 

rib = r X (Jb X B) = (r X f&)- 

IV. THE TIME-DEPENDENT CASE 


Consider an ideal (point) electric/magnetic dipole— 
its total charge is zero, but it carries an electric dipole 
moment p(t) and a magnetic dipole moment m(t). Its 
position (r') is fixed, but its dipole moments vary in mag¬ 
nitude and/or direction. It produces scalar and vector 
potentiali^ 


1 ^ 

V(y,t) = - -2 • 

47reo ^ ^ 

A(r,t) = 

^ ’ 47r| ^ 


P(tr 


-P(^r) 


(42) 


m(i^) H-m(ir) 

c 


(43) 


where the dots denote time derivatives, and the sources 
are evaluated at the retarded time 




(44) 


3-dep' 

larization and magnetization are thereforeI21 


The potentials of an object with time-dependent po- 

d^r'. 


V = 


A=^ 


1 

47reo 

do f 

47r ./ 


n. 

P(r',tr) 


P{r',tr) H- P{r', tr) 

c 


(45) 


M(r',tr) H-M(^^^^) 

c 


As always, we use Eq. and integrate by parts: 
1 


>dh'. 


V = 


ineo 


V' 

1 


— (P+—P 

^ V C 


d^r', 


V' ■ 


n- ■ 

P-f P 

C 


d^v' 


(46) 


where V • P denotes the divergence with respect to the 
the first argument (the explicit r') only. Now, from 


(48) 


(40) 

Eq.|44) 

v'u = --V'^ , 

c 

(41) 

and 



= — x')"^ + {y — y'Y + {z — z'Y 


= -n, . 


(49) 


So 


V' • 


■ 


~ i ■ 

p + —p 

= 

V P+— P 

c 


c 


n, ■ n- . ■ 

-. p -I -V P 

c c 

= V P+—V' P. 


(50) 


Meanwhile, the first term in Eq. can be converted to 
a surface integral: 


1 


P • n j , , 1 


V = d) ^^da'+- I V'-Pd^r' 

47 reo \ Js ^ c dv 


ETAl.V.p 


d‘r’}. (51) 


The two P terms cancel, and we are left with 


E(r,t) = 


47reo 


Pb(r'iir-) .3 , , / cr6(r',t^) , , 


-d"r' + 


Uv 


-da 


(52) 

The bound charges are unchanged (Eq. , though they 
are evaluated, now, at the appropriate retarded times. 
Turning to the vector potential (Eq. 1^ 


A= / —d^r'- [ V X 


47r 


— (M-f —M 
^ V c 


d^r' 


— V' X ( M-h —M 
^ V c 


dV . 


(53) 


Proceeding as before, 
V' X 


M-h —M 
c 


= V xM + —V' X M, (54) 
c 


and 


-/’'>'(f'*' 




-d^r' 


Note that V' acts not only on the explicit r' dependence 
in P(r',tr), but also the implicit r' dependence in t^. 
Thus 


V' P = V P + P V't^, 


(47) 


do I /" Jb(r , tr) + Jp(r , tr) ^3^1 

47r I dy ^ 


( 55 ) 
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where 

5P 

j, = a ^ <“> 

Again, the bound currents are unchanged (though they 
must now be evaluated at the retarded ti mes) , but they 
are joined by the polarization current (Eg. |56[ ). 

Next we calculate the force on the polar¬ 
ized/magnetized object. To begin with, we need 
the force on point dipoles (p(t) and m(t)), in the 
presence of time-dependent fields. The Lorentz force law 
says 

F = [[pE + iJx B)] d^r. (57) 


The charge and current densities for point dipoles (at the 
origin) ar^Sl 

p(r,f) =-(p • V)(5^(r), (58) 

J(r, f) = p(5^(r) — (m X V) ^^(r), (59) 


F = I { - [(P • V)53(r)]E -f [p53(r)] x B 

— [(m X V)(5^(r)] x B|d^r (60) 

= (p • V)E -I- (p X B) -I- m X (V X B) -I- (m • V)B, 


where E and B are evaluated at the location of the 
dipoles. Except for the addition of the p term, the force 
on time-dependent dipoles in time-dependent fields is un¬ 
changed from the static case (Eqs.j^and 22). 

The total force on a chunk of polarized/magnetized 
material is thus (integrating by parts as in Eg. and 
going through steps similar to those leading from Eg. 
to Eg. 28) 


F = 


(P • V)E -h (P X B) -f M X (V X B) 


(M-V)B 


d-^r 


[ [ph E -|- (Jfe -|- Jp) X B] d^r 
Jv 

-f ^ [cijE -|- (Kb X B)] da. 


(61) 


(62) 


This is precisely the force acting on the bound 
charges/currents and the polarization current. As al¬ 
ways, the bound quantities get the total force right. But 
the force density suggested by Eg. 

fb = (-V • P)E-h (P X B)-k (V X M) X B, (63) 


is not at all the same as the actual force 


density (Eg. 61) 


f = (P-V)E-|-(PxB)-|-Mx (VxB) + (M-V)B. (64) 


The torque on a (time-dependent) electric/magnetic 
dipole iJi^ 

N = (p X E) -I- (m X B) -I- r X (65) 

[(p • V)E -f (p X B) -I- m X (V X B) -I- (m • V)B]. 


The total torque on a piece of polarized material is there¬ 
fore 

N = y |(P X E) -h (M X B) 

-fr X [(P • V)E-h (P X B) (66) 

+ M X (V X B) -h (M • V)B] } d^r, 

or, integrating by parts as before (Eqs. 15 and 39): 

N = y r X [ffbE -I- (Kb X B)] da 

I r X [pbE -)- (Jb -f Jp) X B] d^v. (67) 
Jv 


Equation 66 says the torque density is 

n = (P X E) -I- (M X B) -I- r X f 


( 68 ) 


(where f is given by Eg. 64) but Eg. 67 suggests a differ¬ 
ent torque density 


nb = r X fb 


(69) 


(where fb is given by Eg. 63). 


V. THE EINSTEIN-LAUB FORMULA 

The fundamental force law in classical electrodynamics 
is 


F = g[E-k(vxB)], or f = pE-k (J x B) (70) 

(known universally as the “Lorentz force law”). If you 
separate the charge and current into free and bound 
parts, 

p = p/ -f Pb = P/ - V • P, 

J = J^+Jb+Jp=J/ + (VxM) + ^, (71) 

and substitute this in, you get Eg. (including now any 
free charge/current terms): 

fu = p/E -k (J/ X B) 

- (V • P)E-k (V X M) X B-k (P X B). (72) 

In the optics community Eq.j^is sometimes itself called 
the “Lorentz force law” (that’s why we use the subscript 
L).l^This terminology is misleading. As we have seen, the 
substitution (Eg. is incorrect when calculating force 
and torque densities, though it does (when combined, 
of course, with the appropriate surface terms) yield the 
right total force and torque on an object. By contrast, 
Eg. treats the material as a collection of electric and 
magnetic dipoles, not as a distribution of bound charges 
and currents: 

f = P/E + (J/ xB) + (P-V)E 

-f(PxB)-f (M-V)B-kMx (VxB). (73) 
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The fact that their integrals are equal suggests that Tl 
and f differ by a total derivative. Indeed, 

f - 4 = (P • V)E + (V • P)E + (M • V)B 

+ M X (V X B) - (V X M) X B 

= V(M • B) + [(P • V)E + (V • P)E] 

-[(B-V)M + (V-B)M]. (74) 

Now 

[(P • V)E + (V • P)E], = + {VjPj)E, 

= V,(P,E,), (75) 

(and similarly for M and B), so 

(f - fn)* = y^iM,Bj) + V, [PjE, - M,Bj]. (76) 

There is a final twist to the story. By “force” we mean, 
of course, the rate of change of momentum. But in spe¬ 
cial relativity the momentum of a system consists of two 
parts: “overt” momentum associated with motion of the 
center-of-energy, and “hidden” momentum,Ii^ associated 
with internally moving parts but not reflected in motion 
of the system as a whole. Thus 


Ph- 


(77) 


If we are only interested in the overt motion, we might 
introduce an “overt” force, 


_ dpo dph 


(78) 


Now, the hidden momentum of a magnetic dipole in 
an electric field i^ 


Ph = 4(m X E), 


(79) 


so the overt force density on magnetized material is 

1 9(M X E) 


f„ = f- 


dt 


(80) 


Thus 


fo = P/E + (j/ X B) + (P • V)E + (P X B) 
-h M X (V X B) -h (M- V)B 

-4(MxE)-4(MxE). 

This is almost the “Einstein-Laub” force density,li^ 


(81) 


fEL = P/E + [j/X (poH)] + (P • V)E + P X (poH) 
-h (M • V)poH - 4 m X E 
= p/E + (J/ X B) -b (P • V)E + (P X B) 

-b (M • V)B - po f(J/ X M) -b (P X M) 


(M- V)M-beo(M X E) 


(82) 


In fact, using 

V X B = PqJ 4— 
c"' 

= Po^J/ + P + V X M -b cqE^ , (83) 

we get 


fEL = fo-yV(M2). (84) 


Since the “extra” term (—(po/2)V(M^)) is a pure gra¬ 
dient, it will not affect the total force on an object-but 
it does, of course, change the force density^ 

The same considerations apply to torque: the total 
angular momentum consists of two parts, 

L = L„ + L;,. (85) 


The overt torque is 


dLo dL dL/i 
dt dt dt ’ 


( 86 ) 


where 


L;, = r X = 4r X (m X E) (87) 

is the hidden angular momentum of the magnetic 
dipole. Thus the overt torque density on polariz¬ 
able/magnetizable material is (cf. Eq. [6^ 


rio = (P X E) -b (M X B) -b r X fo. (88) 


Meanwhile the Einstein-Laub torque density is 


hel = (P X E) -b (M X B) -b r X f] 


EL 


= no- X (VM^). 


(89) 


Notice that rio and hel yield the same total (overt) 
torque on an object, though they describe rather different 
torque densities. 

In recent years some authorJ^ have advocated the 
Einstein-Laub force law (Eq. [8^ , as a replacement for 
what they call the “Lorentz” law (Eq. 72). We agree 


that the latter is defective, but proponents of the former 
should be aware that they are only talking about the 
“overt” part of the force density, and including an extra 
term (Eq. 841 of dubious provenance. 


VI. CONCLUSION 

So, what is the use of bound charge (and bound current 
and polarization current)? When is the substitution 

p = p/ -b p& = p/ - V • P, 

i9P 

J = J/ + Jb + Jp = J/ + (V X M) -b “y I 
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(Eq. 711 legitimate? Answer: it’s fine for calculating po¬ 
tentials and fields, and hence for use in Maxwell’s equa¬ 
tions. It’s OK when you are interested in total forces and 
torques. But it does not yield the right force and torque 
densities —it distributes the force (over the object) incor¬ 
rectly, even though it gets the total right. There is noth¬ 
ing wrong with the Lorentz force law (Eq. 70) itself.^ 
The problem, rather, is that the substitution Eq. [71] does 


not take proper account of the “internal” surface bound 
charge and current. 
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^ Surface bound charges (and currents) are implicitly in¬ 
cluded in the volume terms, as delta-function discontinu¬ 
ities at the boundary, but it is generally safer to handle 
them as separate contributions. 

^ See, for example, D. J. Griffiths, Introduction to Electrody¬ 
namics, 4th ed. (Pearson, Upper Saddle River, NJ, 2013), 
Sections 4.2.1 and 4.2.2. 

® In an ionic lattice the definition of P is ambiguous, and 
so is the distinction between free and bound charge. See 
E. M. Purcell and D. J. Morin, Electricity and Magnetism, 
3rd ed. (Cambridge University Press, Cambridge, UK, 
2013), Section 10.14. In this paper we shall assume the 
individual dipoles are unambiguously identifiable. 

^ There is a vast and contentious literature on this subject. 
See, for example, P. Penfield and H. A. Haus, Electrody¬ 
namics of Moving Media (MIT Press, Cambridge, MA, 
1967), S. R. de Groot and L. G. Suttorp, Foundations 
of Electrodynamics (North-Holland Pub. Co., Amsterdam, 
NL, 1972), A. C. Eringen and G. A. Maugin, Electrodynam¬ 
ics of Continua (Springer-Verlag, New York, NY, 1990), 
S. Bobbio, Electrodynamics of Materials (Academic Press, 
San Diego, CA, 2000). 

® Actually, there is some subtlety to this, when applied to 
points inside the material (see Section 4.2.3 in ref. 2, or 
Section 2.6 of Bobbio, ref. 4). 

® See, for example, ref. 2, Sect. 4.1.3. 

^ At a surface charge (or current) the electric (magnetic) 
field is discontinuous. The standard limiting argument says 
the correct field to use in such cases is the average 


E — 2 (^outside -f Ejnside), 


B= |(Bo 


+ Bi, 


See, for example. Section 1.14 of Purcell and Morin, ref. 3. 
We take the surface integral in Eq. 9 as the limit f —>■ 0 of 
the volume integral of pi,E over a thin layer of thickness t 
at the object’s surface, in which the polarization P drops 
continuously from its interior value to zero. 

® This is sometimes called the “Kelvin” force. We will use no 
subscript for densities based directly on the dipoles, and a 
subscript h for densities based on bound charge. 

® M. Mansuripur, “Electric and Magnetic Dipoles in the 
Lorentz and Einstein-Laub Formulations of Classical Elec¬ 
trodynamics,” Proc. SPIE vol. 9370, 93700U~1-15 (2015) 
[arXiv:1503:02111 (2015)]. 

See, for example, ref. 2, Eq. 6.3. This formula is often writ¬ 
ten in the equivalent form F = V(m ■ B), but that could 
be misleading in the present context, since m is a fixed 
quantity, and is not to be differentiated. 


See, for example, ref. 2, Sect. 6.1.2. As in Eq. 12 the first 


term is the torque about the center of the dipole, and sec¬ 
ond is (r X F), where F is the force on the dipole. 
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